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Abstract. The possibility of detecting weak lensing ef- 
fects from deep wide field imaging surveys has opened new 
means of probing the large-scale structure of the Universe 
and measuring cosmological parameters. In this paper we 
present a systematic study of the expected dependence of 
the low order moments of the filtered gravitational local 
convergence on the power spectrum of the density fluc- 
tuations and on the cosmological parameters f2o and A. 
The results show a significant dependence on all these pa- 
rameters. Though we note that this degeneracy could be 
partially raised by considering two populations of sources, 
at different redshifts, computing the third moment is more 
promising since it is expected, in the quasi-linear regime 
and for Gaussian initial conditions, to be only f2 depen- 
dent (with a slight degeneracy with A) when it is correctly 
expressed in terms of the second moment. 

More precisely we show that the variance of the con- 
vergence varies approximately as P(k) f2g' 5 z\~ 5 , whereas 
the skcwness varies as 8 z- 1 - 35 , where P(jfe) is the pro- 
jected power spectrum and z s the redshift of the sources. 
Thus, used jointly they can provide both P(k) and f2o- 
However, the dependence on the redshift of the sources is 
large and could be a major concern for a practical imple- 
mentation. 

We have estimated the errors expected for these pa- 
rameters in a realistic scenario and sketched what would 
be the observational requirements for doing such measure- 
ments. A more detailed study of an observational strategy 
is left for a second paper. 

Key words: Cosmology: Dark Matter, Large-Scale Struc- 
tures, Gravitational Lensing 
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1. Introduction 

The detection of weak gravitational lensing effects (weak 
shear and magnification bias) at large distance from clus- 
ter centers is now well established (see Fort & Mellier 1994 
and Narayan & Bartelmann 1996 for reviews). With the 
present development of wide field CCD detectors mounted 
on the best telescopes it becomes possible to start the 
investigation of large scale structures from lensing effect 
observations in an angular area of significant size. The 
shear and the magnification are expected to be imposed 
on the background galaxies by the density fluctuations 
along the line-of-sight. Compared to usual determinations 
of the density fluctuations in the local universe with galaxy 
counts (in fact, the light distribution), this approach is ex- 
tremely attractive since it is directly sensitive to the global 
mass fluctuations regardless of any biases associated with 
the light distribution. The detection of weak lensing has 
thus naturally been suggested as a possible means to mea- 
sure the power spectrum of the large scale density fluctu- 
ations (Blandford et al. 1991). 

Theoretical studies done by Blandford et al. (1991), 
Miralda-Escude (1991a), Kaiser (1992) and more recently 
by Villumsen (1996) proposed various approaches and dis- 
cussed the feasibility of observational strategies for some 
specific cosmological scenarios. In all cases, they concluded 
that weak lensing induced by large scale structures is de- 
tectable with the present techniques. They claimed that 
the projected power spectrum should be recovered pro- 
vided some crucial observational issues, namely the cor- 
rection of image degradation, are solved. However, even if 
the observational aspects are certainly the main difficul- 
ties in the future, different theoretical problems have not 
been addressed yet. In particular, the dependence of the 
physical quantities (such as the convergence, see hereafter) 
on the cosmological parameters have not been discussed in 
complete detail. Moreover the errors associated with those 
quantities for realistic scenarios of large scale structure 
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formation. This paper thus is a theoretical study in which 
we explore the dependence of a priori physical quantities 
on the cosmological parameters. We will focus our analy- 
sis on the second and third moment of the local filtered 
convergence at large angular scale, i.e. in a regime where 
Perturbation Theory results are expected to be valid. We 
take here advantage of many results that have been ob- 
tained in the last few years in this regime, showing that 
the Perturbation Theory calculations give extremely ac- 
curate results for the behavior of the high order moments 
of the cosmic density at large scale (many papers should 
be quoted, see Juszkiewicz & Bouchet 1995 or Bernardeau 
1996b for short review papers on the subject). 

Our aim is to present quantitative predictions on the 
behavior of the moments of the one-point probability dis- 
tribution function of the local convergence as a function 
of the matter power spectrum, the cosmological param- 
eters fi and A. As long as one is interested in its sec- 
ond moment in the linear regime, results can easily be 
extrapolated to other quantities such as the magnification 
or the distortion (see section 3 for more details on these 
quantities and their relationships). For the third moment 
however, the choice of the local convergence (or any other 
scalar quantity) is crucial. The distortion is intrinsically 
irrelevant since its third moment, for obvious symmetry 
reasons, should vanish. The first non-trivial moment would 
then be its kurtosis, the connected part of the fourth mo- 
ment, but this is then more intricate to calculate, and also 
would be more difficult to measure. Other scalar quantities 
could have been considered, such as the local magnifica- 
tion, but it is not (at the second order) directly propor- 
tional to the local density, as it is the case for the conver- 
gence. This is why we have preferred this latter quantity. 

The paper is structured as follows. In §2 we present the 
basic equations, valid for any cosmological model, from 
which we derive the quantities of interest. In §3 we ex- 
plore the dependence of the variance of the convergence 
on the cosmological parameters. In §4 we consider third 
order moments. In the last section we present a study of 
the expected errors due to the use of a finite sample that 
are expected to affect the measurements of the second and 
third moment. The errors have been quantified for a realis- 
tic scenario (APM like power spectrum in an Einstein-de 
Sitter Universe). We eventually propose, in view of the 
previous results, an adapted strategy to do such obser- 
vations. Throughout the paper we adopt the convention 
Hq = c = 1, the physical distances have thus to be multi- 
plied by 3000 for H = 100 /cm/s/Mpc to be expressed in 
Mpc. 

2. The physical model for the galaxy and mass dis- 
tributions 

The gravitational lensing effects makes intervene both the 
rcdshift distributions of the lenses and the sources. In the 
last section we discuss in more detail the knowledge we 



have or need to have on the sources, but at this stage we 
simply describe the redshift distribution of these objects 
by n(z)dz with the normalization, 



Azn(z) = 1 . 



(1) 



To illustrate the results we will either assume that the 
sources are all at the same redshift or are given by a broad 
distribution. In the latter case we adopt a reasonnable an- 
alytic model which reproduces the redshift distributions 
observed in the faint redshift surveys and those expected 
from models of galaxy evolution (Chariot & Fall in prepa- 
ration), 



n(z)dz = 



r(9/8) 



exp [-(z/zsf 



Az 



(2) 



in which it is assumed that the average redshift of the 
sources is z s and that the width of the distribution is about 
0.4 z s (sec fig. 1, thin lines in the upper panel). 

We are interested in the lensing effects induced by 
the large-scale density fluctuations that are assumed to 
have emerged from Gaussian initial conditions. In the lin- 
ear regime we thus assume that the local density can be 
written 1 , 



<5(x, z) 



-I 



d 3 k 

(2 tt) 3 / 2 



D+(z) <5(k) e 



ikx 



(3) 



where D + (z) is the time dependence of the linear grow- 
ing mode. Note that the function D + (z) depends on the 
cosmological parameters as well. It is proportional to the 
expansion factor only for an Einstein-de Sitter universe. 
The Fourier transforms of the initial local density contrast, 
<5(k), are assumed to form a set of Gaussian variables. In 
such a case the power spectrum P(k) entirely determines 
the statistical properties of the initial density field. It is 
defined by 



(<J(k)<y(k')> =<5 Dirac (k + k') P(k). 



(4) 



From the observed density fluctuation in the APM galaxy 
survey it is possible to construct a realistic power spec- 
trum P(k) for both its magnitude and its shape. Using 
the results of Baugh & Gaztahaga (1996) we can use the 
shape (we have set c = H = 1), 



P(k) = 1.2 lO- 



Il + (k/k c y 



i3/2- 



with 



k, 



L50, 



(5) 



(6) 



1 Note that in case of a non flat Universe, the product k • x 
should be understood as a covariant product, k ■ x = \k r + 
©o k± where Do and x are respectively the angular distance 
and the distance along the line-of-sight (§3.3). 
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that reproduces the observed density fluctuations in the 
linear regime 2 . Baugh & Gaztanaga compared the ob- 
served fluctuations in the APM galaxy survey with the 
prediction of this power spectrum and found a good agree- 
ment in the linear regime. Results of numerical simulations 
also show that the nonlinear evolution of the density fluc- 
tuations induced with such a spectrum are in good agree- 
ment with the observed full shape of the angular two-point 
correlation function. 

In this paper, and contrary to previous studies, we do 
not make any assumption on both the density of the Uni- 
verse fio and on the cosmological constant, A. 

3. The Weak-Lensing Equations 

3.1. Description of the Deformation in the Image Plane 

The aim of this subsection is to present the mathematical 
objects that describe the local gravitational deformation. 
In particular we recall the relationships between the local 
convergence k, the statistical properties of which we will 
investigate, and the quantities that are more directly ob- 
servable. We consider two neighboring geodesies C and CJ 
which are enclosed by a light bundle, and we call 0, the 
two dimensional position angle between C and CJ at the 
observer position. Let us define a. as the angular diame- 
ter distance between C and C at a given rcdshift z. We 
assume that the small angle approximation is valid, and 
we have, 



ai(z) = a T>ij(z) 6j, 



(7) 



where in the case of an homogeneous Universe, T>ij = 
Dq(z)S^ , a being the expansion factor, Sfj the Kronecker 
symbol, and T>o(z) the physical distance, (see definition 
hereafter). In an inhomogencous Universe T>ij describes 
the deformation of the light bundle produced by the mat- 
ter distribution. The angular distance a corresponds to 
an angular position on the sky £ such that, 



a 

aV 



(8) 



The shapes of the light bundles on the image plane are 
described by the amplification matrix A, the Jacobian of 
the transformation from a virtual screen located at z (the 
source plane), to the image plane as seen by the observer. 
Its inverse, -4 _1 , is actually more closely related to ob- 
servable quantities, 



A 1 i0 {a,z) = 



(9) 



2 We are aware that the adopted power spectrum does not 
have the expected behaviour, P(k) ~ fc~ 3 , when k — ► oo. But 
this should not be a problem since we focus our calculations 
on the linear or quasilinear regime. The large k behaviour of 
P(k) is therefore not relevant for the purpose of this work. 



The amplification matrix can be expressed in terms of the 
geometrical deformation of the light bundle, the conver- 
gence k and the shear components 71 and 72 (Schneider 
et al. 1992), 



A~ 1 (Z) = 



1 - K - 71 

-72 



-72 
1 - K + 71 



(10) 



In particular the local convergence k(£) is given by the 
trace of this matrix, through the equation, 



«(£) = l-tr [A-\M)] /2. 



(11) 



The intensity of the shear 7 = 71 + i-f 2 is given by the 
ratio of the eigenvalues of the matrix A . The relations 
between the physical quantities (k, 7) and the observable 
quantities, the magnification (=amplification) u. and the 
distortion 5 of the images, are given by, 



1 



1 



// 



\A~ 



(1 — k) 2 — 7 1 



; 6 = 



2g 



i + l<?l 2 ' 



(12) 



where g = 7/(1 — n). Note that in the weak lensing regime 
we have the simple relations, 



H = 1 + 2k ; 8 = 27, 



(13) 



In practice, the magnification is measured through the 
local density of objects (Broadhurst et al. 1995, Broad- 
hurst 1995, Fort et al. 1996) or through the image size of 
galaxies in different surface brightness slices (Bartclmann 
& Narayan 1995), whereas the distortion is measured di- 
rectly from the shape of the background galaxies (Bonnet 
al. 1994, Fahlman et al. 1994, Smail et al. 1994, Squires et 
al. 1996a, b). One should be careful when order moments of 
the local convergence higher than the variance are consid- 
ered, since non-linear couplings exist in general between 
u. and 6, and Eqs. (13) cannot be used. It is possible to 
get the convergence from Eqs. (12) by the measurement 
of both the magnification and the distortion, or by using 
only the distortion and the relation introduced by Kaiser 
(1995, see also Seitz & Schneider 1996) which is valid even 
beyond the linear regime, 



Vlog(l-K) = 



1+51 

92 



92 

l-9i 



( ggi dg 2 

d0! d9 2 
-dgi , dg 2 



.(14) 



V dd 2 



+ 



001 



In the following we will thus assume that the local con- 
vergence is accessible to the observations in a given sam- 
ple, and we will focus our investigations on the statistical 
properties of this quantity. 

3.2. The Source Equation for the Deformation Matrix 

We perform the calculations assuming that the Born 
approximation is valid, i.e. that the deformation of a 
light bundle can be calculated along the unperturbed 
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geodesies. This assumption will be discussed in section 
5.7. In a Friedmann-Robertson- Walker (FRW) Universe, 
such a geodesic may be parameterized by the time variable 
A, defined by 



dA = - 



da 

ma) 1 



A(a = a ) = 0, 



(15) 



where H(a) is the (time dependent) Hubble constant 
(in units of Ho). Using the Born approximation and the 
geodesic deviation equation, the evolution equation for the 
angular diameter distance a in a lumpy Universe, along 
the unperturbed ray, as a function of A is (Seitz 1993), 



a(\)=K(\)a(\), 



(16) 



where 1Z is the 2x2 symmetric tidal matrix, and depends 
on the second derivatives of the Newtonian gravitational 
potential (Sachs 1961, Seitz 1993), 



K = 
R = 
F = 



R - TZe(F) Jm(F) 
lm(F) R + Ue(F) 

~~^Ra()k k@ , 



(17) 

(18) 
(19) 



where R a p is the Ricci tensor, C a py5 is the Weyl tensor, k a 
is the wave- vector of the light ray and e a is a complex null 
vector propagated along the geodesic such that e a k a = 0. 

Let ip be the 3-dimcnsional Newtonian gravitational 
potential. The comoving Poisson equation is, 



Aip = AttG p 5, 



(20) 



where 6 = (p — p)/p is the density contrast, A is the 3- 
dimensional Laplacian operator, and the derivatives are 
done with respect to the proper distance. Note that Eq. 
(20) docs not depend on the cosmological constant A. 
Then, it is straightforward to show that, 



R = -4nGa 2 p-a 2 (d n <P + d 22 p) 
F = -a~ 2 (dntp - d 22 tp + 2idi 2 tp). 



(21) 
(22) 



From Eqs. (7) and (16) the central equation driving the 
matrix T> is then, 



d 2 [ a (z)P„(g,z)] 
dA 2 



= a(z) 1Z ik (£,z) V k j{£,z), 



with the boundary conditions, 
,'dVi 

(»tf)z=o = 



dA 



1. 



(23) 



(24) 



The first equation expresses the focusing condition at the 
observer, and the second the Euclidean properties of the 
space at small redshift. IZij is a symmetric matrix that 



can be written in terms of ip, and p using the expressions 
of R and F, 

^--^.-(J ;)-*-(£ £)• (25) 

Since 4irGp — 3/2 flo/a 3 , for a completely homogeneous 
universe, the matrix TZij is given by, 



and we are left with the equation for the distance 
d 2 [a(z) V (z)} 



dA 2 



= -3/2 (1 + z) 4 n V (z). 



(26) 



(27) 



We know that the solution of this differential equation is 
given by 



V (z) 



with 



Vl-fto-A 



sinh 



Vl-fio-A / 
Jo 



dz' 



E(z') 



. (28) 



E[z] = y/A + (1 + z) 2 (l-n -A) + (1 + z) 3 fio. (29) 

This equation has a known analytical solution only when 
the cosmological constant A is zero. Otherwise it has to 
be integrated numerically. 

3.3. The Expression of the Local Convergence 

In this section we restrict ourself to the linear regime. We 
then assume that the difference between the local value 
of lZij(£,z) and its value for a homogeneous universe is 
small, so that, 



n ij {i,z)-nf{z) 



(30) 
(31) 



where we have introduced an effective potential <f> such 
that, 



\^ {1 \tz) 



where <5mass(£ ; z) is the local mass over-density in the lin- 
ear regime. We assume the small angle deviation approx- 
imation is valid. In that case, the plane-parallel approxi- 
mation works, which means that only waves perpendicular 
to the line-of-sight contribute to lensing, and consequently 
we can neglect the second order derivative along the un- 
perturbed ray in the Laplacian operator. Moreover, the 
difference T>ij{^, z) — T>o(z) Sfij is a small quantity, which 
can be expanded with respect to the initial density. We 
then define T>\j\^, z) as the part of this difference which 
is linear in this density field (terms of higher order will be 



8^ 

mass 



(€.*) 



A-KGp'' 



(32) 
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considered in section 5). The differential equation for 
can then be derived from (23) and it reads, 



-a(z) U^{z) V%&z) = 



dA 2 

-|fi (l + zf V (z) 0<J>(€,z) 



with 

P>V) X=0 = , 



dV, 



(i)' 



0. 



(33) 



(34) 



2 = 



To solve this differential equation it is more convenient to 
write it with the variable z. The differential equation then 
reads, 



dz 2 

1 1 dE(z) 



1 dE{z) dVW(t, 



1 + z E(z) dz : ' ' 

3 fi (1 + z) 



+ 



E(z) dz dz 

3 n (i + z) 

2 E 2 {z) 11 



Y" U E 2 (z)~' Vo ^ Z)(t, y {tz) 



(36) 



The homogeneous differential equation associated with 
this equation has two known solutions, the distance func- 
tions T> (z) and Uo(z) are given respectively by Eq. (28) 
and, 



Uo{z) = 



-v/l-fio-A 



cosh 



. (37) 



from which the solution of the differential equation (33) 
can be written. Indeed the solution of the inhomogeneous 
differential equation reads 



; (i+zr Powg(z') 

dz TT-— — — x 



E 2 (z') 

U (z)V (z')-U (z')-D (z) 
U^z>)V {z')-Uo{z>)V> Q {z>y 



(38) 



which, after elementary mathematical transformations 
simplifies in, 



(39) 



sinh 



Vl-fio-A^ ^] (1 + z>) Vo(z') 



It follows from Eq. (11) that the local convergence can be 
written for a source at redshift z (see also Seljak 1996), 



« (1) «) = ~< 

X>o(z',z) Pp(z') 
T> (z) 



! / d x (z')x 
Jo 



where T>o(z', z) is the distance between the redshifts z' and 
z. We have also introduced the new variable x{ z )i very 
useful for such calculations, which is the physical distance 
along the line-of-sight, 



d X (z) = 



dz 

wr 



(41) 



It coincides with T> (z) only for a flat geometry f^o+A = 1. 
Eq. (40) expresses the fact that the local convergence 
A/i = 2kW is given by the superposition of the conver- 
gence induced by each lens between the observer and the 
source. For the distribution of sources n(z) we thus have 



.(i) 



(€) 



dz 



E(z') 



-w(z')5£l ss (£,z'), 



with 

w(z') = 



f> (1 + T> (z') / 

J z' 



dz n(z) 
D (z) 



(42) 



(43) 



y/l-ila-A 



sinh 



Jz> 



dz" 



E(z") 



The previous two equations are the basic ones that relate 
the local convergence to the linear cosmic density along 
the line-of-sight for any cosmological model. 

3.4- The Efficiency Function 

We call w(z) in Eq. (43) the efficiency function. It de- 
scribes the efficiency with which a lens at a given redshift 
z located along the line-of-sight will affect the local conver- 
gence. This efficiency function obviously depends on the 
redshift of the sources. It is maximum at half the distance 
between the source and the observer, and vanishes at both 
ends. In fig. 1, we present the shape of the efficiency func- 
tion in different case. In la, this is for an Einstein-de Sit- 
ter Universe for two different populations of sources. The 
thin lines show the shape of the distribution functions of 
the redshifts of the sources that are assumed to be either 
centered on z = 1 (solid line) or z — 2 (dashed line). 
The typical redshift of the lens is about 0.4 to 0.5, with 
a broader distribution when the redshift of the sources is 
larger. 

In the lower panel we show the dependence of the effi- 
ciency function on the cosmological parameters. Here the 
sources are simply assumed to all lie at z = 1.5. We can 
see that the shape of the efficiency function is not re- 
ally affected. The amplitude is however changed, and it 
is roughly proportional to fio- The dependence on A, al- 
though not totally negligible, is considerably weaker than 
that on flo. 



(40) 

4. The Variance of the Smoothed Convergence 

We are interested in the statistical properties of the local 
convergence at a given scale. We are more particularly 
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Fig. 1. Shapes of the efficiency functions in different cases. In 
the upper panel the selection functions n(z) are given by the 
thin lines (corresponding to Eq. [2] for z a = 1 and z s —2) and 
we assume an Einstein-de Sitter Universe. In the lower panel 
the sources are assumed to be all located at z s — 1.5 and the 
cosmological parameters are varied, thick solid line for £1 = 1, 
A = 0, thin solid line for £1 — 0.3 A = and thin dashed line 
for £l = 0.3 and A = 0.7. 



interested in the variance at a scale corresponding to the 
linear or the quasi-linear regime where the interpretation 
in terms of the cosmological models is easier. For the usual 
cases, sources at a redshift of about unity, it corresponds 
to angular scales of about 30', or more. 



4-.1. The Convolution with the Window Function 

Throughout this paper we assume that the geometrical 
filtering is done with a top-hat window function of angu- 
lar scale 0o- In the expression (42) the smoothing applies 
to the exp(ik • x) factor and it gives in the small angle 
approximation (see Bernardeau 1995), 



—r / dO I d(f> exp(ik • x) = 
1" #o Jo Jo 

W 2D [k ± V (z)9 } exp ifc '*W. 



(44) 



where k r is the component of k along the line-of-sight, kj_ 
its component in the transverse directions and W 2 d is the 
Fourier transform of the 2D top-hat window function, 



W 2D (k) 



' k ■ 



(45) 



( Ji is the Bessel function of the first kind). Note that the 
expression (44) involves both the distance T>q(z) and the 
variable x(z). These two quantities are identical only for 
cosmological models with zero curvature. 

4-2. The Variance in the Linear Regime 

Using the definition (3) of the power spectrum and the 
relations (44, 42) we get 



dzi 



dz 2 



w(z 1 ) w(z 2 ) D+(zi) D + (z 2 ) 



E( Zl )J E{z 2 ) 
d 3 k 



(46) 



(27T) 



P(k) x 



W$ D [h± D (z) 6 ] exp [ik r (x(^i)-x(^))]. 

where D + (z) is the growing mode of the density contrast 
(cf. Eq. [3]). 

We use again the assumption that the deflecting angle 
and the smoothing angle are small. It implies that 



P(fc) w P(fcx). 



(47) 



Then the integral over k r yields a Dirac delta function in 
x(zi) — x{z 2 ) leading to the expression, 



f 

Jo 



dz 



w{z) 2 D\{z) x 



E(z) 

k± dk ± P(k±) W2 D [k ± D (z)$ ]. 



(48) 



In the following section we explore the properties of this 
function with different hypotheses for the power spectrum, 
the redshift distribution of the sources, and the cosmolog- 
ical parameters. 

4-3. Realistic Results for an Einstein-de Sitter Universe 

For an Einstein-de Sitter Universe, the distance takes a 
simple expression as a function of the redshift, 



V (z) = 2-2/VTTz' 
and we have, 

w(z') = %n (l + z')V (z') x 

Cdzn(z) [l-V Q {z')/V {z)\. 



(49) 



(50) 



We also know that the linear growth rate of the fluc- 
tuation is proportional to the expansion factor so that 



D + (z) = 1/(1 + z). 



(51) 
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Fig. 2. The expected variance of the convergence as a function 
of the smoothing angle (in degrees) for z 3 = 1 (solid line) and 
z s — 2 (dashed line) in (2) and assuming an Einstein-de Sitter 
Universe. 

Using the power spectrum given in (5) and the redshift 
distribution (2) for z s = 1 and z s = 2 we have calcu- 
lated the angular dependence of the variance of the local 
smoothed convergence. The results are presented in Fig. 
2. We can see that the typical amplitude of the conver- 
gence is of the order of a few percent. It starts to bend at 
an angular scale of about 5 to 10 degrees corresponding to 
the physical scale where the power spectrum bends down. 
This scale will be of crucial importance for the evaluation 
of the cosmic errors on the measured moments. 

We can also notice that the results have a significant 
dependence on the adopted redshift distribution of the 
sources. In the next subsection we discuss in more detail 
this point. 

4-4- Dependence with the redshift of the Sources 

We explore here the dependence of the second moment 
with the redshift of the sources assuming that they are all 
at a given redshift z s . We also assume that the power spec- 
trum is given by a power law spectrum. For an Einstein-de 
Sitter Universe the result reads, 



<K 2 (0 O )) oc f 

Jo 



T> (z e ) 



dVV 2 {l-V/V {z s )) 2 V- 



©0 (*-) 



-2-n 



(l-n)(2-n)(3-n) 



(52) 



where the angular distance T>q is given in the equation 
(27). For sources of redshift of the order of 1 or 2, and for 
a power law index of n w —1.5, we have approximately, 



(« 2 (0 o ))cxz s 1 - 5 . 



(53) 



This result confirms the trend observed in the previ- 
ous subsection. It already indicates that the redshift de- 
pendence of the second moment is rather large. Therefore 
any interpretation of a measurement of such a moment in 
terms of magnitude of the power spectrum would require 



a precise knowledge of the sources that have been used. 
Moreover the magnitude of the second moment depends 
as well on the values of the cosmological parameters fio 
and A. The dependence of the results on those parameters 
is explored in the next subsection. 

4-5. Dependence on the Cosmological Parameters 



L0g[<K 2 (9)>] 




Log (ft) 



Fig. 3. The expected dependence of the magnitude of the vari- 
ance on Qo for z s = 1 (solid line) and z s — 2 (dashed line). 
The overall normalization is arbitrary. 



In this subsection we assume that the power spectrum 
is a power law and that the sources are all at redshift z s = 
1 or z s = 2. In Fig. 3 we then plot the dependence of the 
second moment on fi for A = 0. We have approximately 



for 2, = 1 and 



(k 2 (0 o )) cx ft, 



1.4 



(54) 



(55) 



for z s — 2. The resulting 0,$ dependence is thus weaker 
than what one would naively expect from Eq. (48) (e.g. 
(k 2 (# )) oc fig). The fact that the fi dependence is actu- 
ally weaker is mainly due to the growth factor of the linear 
mode. Indeed at high redshift, for a given normalization 
at z = 0, the lower Q the larger the density fluctuations 
are. 

In Fig. 4 we present a contour plot of the fl and A 
dependence of the second moment. We can see that the 
expected magnitude of the second moment depends essen- 
tially on f2o- It is however degenerate with A when fio is 
large. As expected, the A dependence is more important 
when the redshift of the sources is larger. 
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Fig. 4. Contour plots of the expected dependence of the magnitude of the variance on fio (horizontal axis) and A (vertical axis) 
for z s = 1 (left panel) and z s — 2 (right panel). The lines are iso- values of the variance, normalized to the Einstein-de Sitter 
case and regularly spaced in linear scale. The variance is stronger at the top right corner of the panels. 



4-6. The Variance with two Populations of Sources 



Interestingly when one considers jointly two different pop- 
ulations of sources, here with z s = 1 and z s = 2, the ratio 
of the two variances depends on fio and A (see contour 
plot of Fig. 5) but independent of the normalization of 
the power spectrum. (The dependence on the power law 
index is weak). 

This is thus, a priori, a way to disentangle the different 
parameters. 



5. Higher-Order Moments 



In this section we explore the possibility of considering 
higher order moments of the distribution of the local con- 
vergence. The idea is that the higher order moments are 
sensitive to nonlinear aspects of the dynamics and thus 
could provide independent constraints on the cosmologi- 
cal parameters. 

The nonlinear terms in the expression of the fields can 
be calculated using Perturbation Theory. Similar calcu- 
lations have been done for the density field and for the 
local velocity divergence. Comparisons with numerical re- 
sults have shown that Perturbation Theory calculations 
give excellent quantitative predictions of the behavior of 
the high order moments (see Bernardeau 1996b and refer- 
ences therein). 



5.1. Second Order Perturbation Theory for the Local Con- 
vergence 

The basic assumption of perturbation theory is that the 
local convergence can be expanded in terms of the initial 
density field, 



A2) 



(€,*) 



(56) 



where k^(£, z) is linear in the initial density field, so in 
the variables <5(k) (this is the term calculated previously), 
and k( 2 '(£,z) is quadratic in those variables, etc... For 
calculating the expression of (£, z) we need to use the 
differential equation (23) up to its second order introduc- 
ing the expression of 



with 



~Q (i+^) 5 $he,z) 



8 m 

u mass 



(57) 



(58) 



To find the expression of the local convergence at the sec- 
ond order it is thus necessary to have the expression of 
the local density at the same order. 

The calculation of the second order term for the local 
density has been investigated in many papers. It has now 
been calculated for any cosmological model, and a very 
useful description is given by, 

4 2 ) S8 (x) = D\{z) J d^MfaWk,,) x 

k x .k 2 , 2(k!.k 2 ) 2 ^ 

exp [lx • (ki + k 2 )j, (59) 



kl 



h 2 h 2 
«, 1 «, 2 
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where the <5(k) factors are the Fourier transform of the 
initial density field (cf. [3]). 

The coefficients appearing in this expression given for 
an Einstein-de Sitter universe, are in fact slightly Qq and A 
dependent (Bouchet et al. 1992, Bernardeau 1994). These 
dependences are actually so weak (less than 1%) that they 
can be safely neglected in the subsequent calculations. 

(2) 

The differential equation for V\- (£, z) can be obtained 
from the equation (23), when it is developed up to its 
second order. We then get, 



dA 2 

-|n„ (l + z) 4 ^{tz) V Q (z) + 
+a(z) K&\t,z) X>£>(€,z) 



with 



(V 



(2)- 
ij ■ 



dV. 



(2) 



z=0 



(60) 



(61) 



The expression of the local convergence can be deduced 
from the previous two equations, 



* (2) (€) 



E(z 



- w(z')x 



(62) 



€L(tz') ~ 5 tr [^(tz')V^(tz')] /V (z') 



As it can be observed two terms are contributing to this 
expression. One term is given by the second order density 
field, the other one is a combination of the linear order of 
the local density and the linear order in the local ampli- 
fication. The latter corresponds to the nonlinear coupling 
introduced by two subsequent deflecting lenses. In the fol- 
lowing wc will see that the contribution to the third mo- 
ment is dominated by the first contribution. 

5. 2. The Expression of the Local Skewness 

The principle of such a calculation has been developed 
in previous papers for the density or the velocity field. It 
relies on the assumption that the initial density field is 
Gaussian. As a consequence we have, 



(5(k 1 )...«5(k 2p+1 )) = 
and 



(63) 



(5(k 1 )...«5(k 2p )>= II <<5(k 2j --i)<*(k 2j -)> (64) 

permutations j—1 

where the sum over the permutations is made in such a 
way that all possible pair associations are taken into ac- 
count. 
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Fig. 5. Contour plot of the expected dependence of the mag- 
nitude of the ratio of the two variances on fio (horizontal axis) 
and A (vertical axis). 



Then, applying those properties to the initial density 
field one can calculate the leading term for the expression 
of the third moment of the local convergence. Using the 
expansion (56) we have 

(n 3 (0 o )) = <(^+^+...) 3 ) 

= ((>)) 3 }+3((>))V>> + ... (65) 

The first term of this expression is identically zero in case 
of Gaussian initial conditions. The dominant contribution 
to the skewness is thus given by the next-to-leading term, 
3( (kW) k*- 2 -*), that combines both the expression of the 
linear amplification, and its second order expression. 

To compute the previous expression it is worthwhile to 
know that (Bernardeau 1995) 

J d 2 k! y"d 2 k 2 VF 2D (|k 1 +k 2 |) W 2D (k 1 )W 2D (k 2 ) x 



P(h) P{k 2 ) 



ki • k 2 



2 (kx • k 2 ) 2 



p 1 7 pp 
Jd 2 ^ Jd 2 k 2 P(kl) P{k 2 ) X 

IwiMwiM + ^w 2 D (k 2 )^i(k{ 

Using this property we get, 

6 f°° dz' 



(66) 



<^o)) = 



(27r)2y E(z 



- w{z'f D%{z>) 



J kdkWi D (k) P(k) 



(67) 
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This integral can thus be integrated numerically without 
more technical difficulties than for the second moment 3 . 
At this stage it is crucial to notice that the ratio, 
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W)) 
<« 2 (^o)) 2 



(68) 



is expected to be independent of the normalization of the 
power spectrum. In the following we will explore in more 
detail the dependence of this ratio on the different cosmo- 
logical parameters and physical hypothesis. 

To start with let us consider a very simple case in which 
we assume that we have a power law spectrum, that all 
sources are at the same redshift z s and that we live in an 
Einstein-dc Sitter Universe. In such a case we find that 




S3(z s ) 



36 3, 



40 



(n- 
4n(2n - 1) 



1) (n - 2) {n- 3) 2 x 



smoothing angle 
Fig. 6. The expected parameter S3 as a function of the smooth- 
ing angle (in degrees) for z a = 1 (solid line) and z s = 2 (dashed 
line) in (2), the power spectrum (5) and assuming an Ein- 
stein-de Sitter Universe. 



32 <, 



/ 



\ + z s VTTz~ : 
16n(2n-5)(2n-3)(2n-l) (2 - 2/v/T+zI) 
24 48 



(69) _ 



5535 40 



1 + Z S 



for n 



32768 2 - 



-3/2; 



which gives, 

s 3 {z s ) = w -42 for n = -3/2 and z s = 1. (70) 

This is a quantity a priori accessible to a measurement in 
a reasonably large catalogue. In the last section we will ex- 
plore in more detail the expected magnitude of the errors 
for doing such a measurement. Let us start with a study 
of the dependence of this quantity on the cosmological 
parameters. 

5.3. Numerical Results for a Realistic Case 

We present in Fig. 6 the results obtained for the function 
53(^0) for the power spectrum (5), the source distribution 
(2) with z s = 1 (solid line) and z s =2 (dashed line) in 
case of an Einstein-de Sitter Universe. We can see that 
the function S3 is rather flat, with little variation with the 
smoothing angle. 

5.4- Dependence on the redshift 

From equation (69) we easily can visualize the z s depen- 
dence of S3. This dependence is shown in Fig. 7. We can 

3 Note that the property (66) is valid for an angular top-hat 
window function only, and would not be valid for a Gaussian 
window function for instance. For other types of window func- 
tion the dependence of the skewness on the shape of the power 
spectrum is therefore expected to be more complicated. 




1 1.5 2 

source redshift, z s 

Fig. 7. The expected dependence of S3 on the redshift of 
the sources for an Einstein-de Sitter Universe. The curves 
have been plotted for different power spectrum indices, 
n = —1.3, —1.5, —1.9 with the thick grey line for the in- 
termediate value 

see that the result is almost independent of the index n, 
and we have approximately, 



S3{Zs) 



-42 z\ 



-1.35 



(71) 



As for the variance, we found that the dependence of the 
skewness on the redshift of the sources is also quite large. 

5.5. Dependence on the Cosmological Parameters 

From the general formulation we can also explore the de- 
pendence of the results on the cosmological parameters. 
In Fig. 8 we present the fio dependence of the result for 
A = 0. We can see that the dependence is rather large, 
approximatively as 



s 3 (fi ) w -42 n " ai 



(72) 



for z,Ril. The Slo dependence is actually slightly weaker 
when the sources are at higher redshift (when the sources 
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Fig. 9. Contour plot showing the dependence of S3 on Qo (horizontal axis) and A (vertical axis) for z s = 1 (left panel) and 
z s = 2 (right panel). The contour lines are the logarithmic values of S3 regularly spaced in a logarithmic scale. 



LogC-s^ 




Log (fl) 



Fig. 8. The expected Qo dependence of S3 for z s = 1 (solid 
line) and z s — 2 (dashed line). 



are at very small z we recover the natural expectation that 
S3 varies inversely with fio-) 

In Fig. 9 we then present the contour plot showing the 
dependence of the parameter on Oo and A. We can see 
that S3 is significantly less A dependent. 

5.6. Discussion, Simplified Model 

A this stage it is worth to compare the qualitative re- 
sults obtained here with those obtained for an a priori 
similar case: the statistics in an angular survey. In both 
cases indeed the observed quantities are proportional to 
the cosmic density integrated along the line-of-sight. The 
expected properties of the moments of the one-point PDF 
are however qualitatively different. The reason is that for a 
2D angular survey, the density fluctuations are calculated 



with a selection function implicitly normalized to unity, 
whereas the efficiency function in the case of the weak 
lensing is not normalized (see fig. 1). As a consequence the 
local convergence is the summation of the effects of each 
lens plane, whereas the local measured galaxy over-density 
identifies with the average of all contributions present on 
the line-of-sight. The behavior of these quantities is thus 
dramatically different when the depth of the respective 
catalogues is changed. The variance of the local conver- 
gence increases with the depth (as for a random walk, 
the mean distance from the center increases as the square 
root of the number of steps), but it becomes more and 
more Gaussian, so that S3 decreases. On the other hand 
the variance of the local galaxy density contrast decreases 
with the depth, since the density fluctuations tend to be 
smoothed out. The parameter S3 is on the other hand al- 
most independent of this depth. This is exactly what one 
would obtain by considering respectively the summation 
or the average of a given number of quasi-Gaussian vari- 
ables. 

From these general remarks one can build a very sim- 
ple model for the moments of the local convergence that 
reproduces qualitatively the results discussed in the pre- 
vious sections. Let us define uj{z s ) as the average of the 
efficiency function when the sources are at redshift z s , 



/o E(z) 



Then, from eqs. (48) and (67) we simply expect that 



-1 



S3 OC — 



(73) 

(74) 
(75) 
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which show the dependence of these statistical quantities 
on the efficiency function and the depth of the catalogue. 
These relations are certainly crude approximations, but 
explain most of the features encountered previously, such 
as the dependence of the results on f2 and on the redshift 
of the sources. It also suggests the existence of a combi- 
nation of (k 2 ) and s 3 that would be almost independent 
of the redshift of the sources. This combination (which is 
almost the product of these two quantities, since ZJ(z s ) is 
proportional to the distance) would contain information 
of purely cosmological nature. 

5.7. Systematics due to other Quadratic Couplings 

As it can be seen from Eq. (65) the skewness is induced 
by any quadratic coupling in the observed convergence. In 
particular it implies that the convergence should not be 
measured from the local shear using the linear approxima- 
tion. The method proposed by Kaiser (1995), or a similar 
method, that takes into account the nonlinear relation be- 
tween the local shear and the local convergence has to be 
used. 

Spurious couplings can appear also through: 

— nonlinear coupling between two deflecting lenses; 

— coupling between two deflecting lenses due to the in- 
duced displacement of the light path; 

— coupling between the population of selected sources 
and the local convergence; 

— density fluctuations of sources. 

In the following we subsequently examine the impor- 
tance of these different effects. 

5.7.1. Coupling between two Deflecting Lenses 

This coupling is due to the fact that when the light path 
intercepts two lenses the resulting effect is given by the 
product of the amplification matrices. The resulting con- 
vergence Ktot is given by the trace of this product, but 
coincides with the sum of the convergences induced by 
each lens separately only in the linear regime. In general 
it is indeed given by 



Ktot = «A + KB - K A K B -JA-IB 



(76) 



which contain quadratic terms in the local lens densities 
(the indices A and B correspond to quantities associated 
with the lenses A and B respectively). 

Actually this nonlinear coupling between two lenses 
appears in the neglected term in (62). This term induces 
a corrective expression for the second order convergence 
term, 



L 



q fOO J f 



The other source of coupling is due to the failure of 
the Born approximation (i.e. the densities are calculated 
along the unperturbed light path) when two lenses are 
involved. Indeed the convergence induced by a given lens 
is not due to the local projected density at the observed 
angular position if a foreground lens shifts the apparent 
position of the background galaxies. It implies that in (40) 
the local over-density can be taken at the unperturbed 
position only at the linear order but in general should be 
taken at the position £ + S£ where 5% is the displacement 
induced by the foreground lenses. To get the quadratic 
coupling one can simply write, 



€L(£ + z) « z) + V5« (€, z) ■ st 

This term induces a term similar to (77), 



« (2) 

v corr.(2) 



00 dz' 



E(z') 



-w(z') x 



(78) 



(79) 



o E(z") V (z') 

These two terms obviously induce a corrective term for 
the skewness. It can be easily calculated, and we estimate 
it here for an Einstein-de Sitter universe for sources at a 
fixed redshift z s . In such a case the second order corrective 
term for the convergence reads, 



K (2) tt\-- r° {Zs \ v M z s)- v ' I av"v" 

{V'-V") W{V')^{V'') + M{V')^\V'')] .(80) 



f 

Jo 



The resulting corrective term for S3 is 

4- = 3(41. (k (1) ) 2 )/((« (1) ) 2 > 2 > 
so that 4 



(81) 



corr. 
S 3 



T>o(z s ) 



V 



6 



-in + 2 — 

2 V ! V 



dV / dV" 

Jo 

(V Q (z s ) - V'f (V (z s ) - V") (V - V") {V"Y x 
V'V (z s ) (6 V')- 2 - n (Oo V")- 2 - n / 

/ dV (V'f (V Q (z s ) - V'f (9 V')- 2 - n 

Jo 

w 1.55 for n = -3/2, 



(82) 
(83) 



o E{z") V (z') 



independently of the redshift of the sources. It is com- 
pletely negligible compared to the main contribution (69), 
and note that this effect is expected, at first glance, to be 
mainly independent of ilo- 

It might be also worth to note that the resulting de- 
formation field induced by two lenses is no more potential 

4 we use here a similar property as in (66), see Bernardeau 
(1995). 
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(the amplification matrix is not symmetric). As a result, 
the method proposed by Kaiser (1995) to build the local 
convergence from the local shear is expected to fail at this 
level. Once again, this is not expected to be crucial since 
the effects of this coupling are seen to be small. 

5.7.2. Coupling between the Population of Sources and the 
Convergence 

The origin of this effect is the fact that the population of 
sources on which the shear is measured may change with 
the local magnification. In this case 



,( 2 ) 



with 



(I) / 

JO 



dz' dw(z') 



E(z') 8k 



(84) 



dw(z') 
8k 



*' 1+ ^«f4^ x (85) 



Vl-fto-A 



sinh 



Vl-fio-A 



dn 
2 dz" 
, E(z") 



where 8n{z)/8n is the partial derivative of the number 
density of sources at a given redshift with the local con- 
vergence (therefore the local magnification) for a fixed 
normalization. It can be evaluated a priori from the pop- 
ulation of objects used to do the measurements, and on 
the selection procedure. For instance the mean depth of 
the sources is expected to change with the local magnifi- 
cation. We can build a simple model in which we assume 
that the sources are always at a given redshift, but that its 
corresponding angular distance is a function of the local 
magnification, 



V (z s )=V a (1 + ok). 



It would imply, for an Einstein-de Sitter universe, 

8w{z) 



iVl{z) 



8k 



2 V° 



(86) 



(87) 



In such a case the corrective term for S3 can be calculated 
straightforwardly and is given by, 



9a(l-n)/2, 
11 a for n = —1.5. 



(89) 



It is not necessarily completely negligible, depending on 
the value of a, that is the sensitivity of the selected popu- 
lation of sources with the local magnification (this will be 
explored in more details in the second paper). So, one way 
to reduce the effect of this coupling is to select the galaxies 
used to measure the shear from a surface brightness crite- 
rion (Bartclmann & Narayan 1995). In practice however, 
it is probably impossible to suppress all influence in the 
selection procedure of the total luminosity of the objects. 



5.7.3. Coupling due to Density fluctuations of Sources 

Throughout this study we have always implicitly assumed 
that the sources formed a uniform background population. 
This is actually far from being true, and one naturally ex- 
pects the sources to have cosmic density fluctuations. The 
effects of these fluctuations of the measurements of the lo- 
cal convergences strongly depend on the adopted method. 
For a method based on the estimation of the local mag- 
nification by the depletion of the galaxy number density 
(magnification bias), this can be a very strong effect on 
both moments. For methods in which the local conver- 
gence is built from the measured shear, this effect is weak. 
Simple considerations show that it introduces a significant 
coupling only at the level of the fourth moment. This is 
clue to the fact that the populations of sources and lenses 
can be considered as two independent populations so there 
are no direct couplings between their density fluctuations. 

6. Error Evaluations due to the Cosmic Variance 

In this section we investigate the effects of the cosmic vari- 
ance on the precision of the measurements of the second 
moment and on S3. The origin of this error is the finite size 
of the sample in which the measurement could be done. 
Typically we will assume that the sample size is about 
25 deg 2 , which could be accessible in the coming years for 
the now available technologies. In this section we describe 
the measured local convergence, K moas .(£) as the sum of 
the exact local convergence and a noise e(£), 



.«i) = «(€i) + c(€i)- 



(90) 



The random noise e(£j) components are assumed to be all 
independent of each other, independent of the true values 
of the local convergence and to have the same variance e, 



(91) 



It is reasonable to assume that e is of the order of 10~ 3 (see 
Appendix B). We then introduce the notation X which 
is the (connected) geometrical average of the quantity X 
in the sample. These quantities are also cosmic random 
variables and are thus expected to vary from one sample 
to another. Our goal is then to estimate the ensemble av- 
erages and variances of those geometrical averages for the 
quantities of interest, i.e. to know how good they are as 
estimates of the true cosmic quantities. 

To do such calculations we have to introduce the ex- 
pression of joint moments of the local convergence. Luck- 
ily, this is not too complicated! Indeed in the quasi-linear 
regime one can show that, when the separation between 
two cells is large compared to the smoothing angle, we 
have 



(s"(0 o ,*i) K q (6 ,t 2 )) = 
c Pq <k(0 o ,€i) «(0o,€ 2 )) (« 2 (fl )) P+9 ~ 2 , 



(92) 
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where c pq are finite quantities. This is the application to 
2D statistics of what was developed in 3D by Bernardeau 
(1996a). However contrary to the 3D case, and due to pro- 
jection effects, these coefficients depend on the details of 
the adopted model. They can anyway be calculated nu- 
merically (see Appendix A). Using these results we can 
calculate both the expectation value of a measured quan- 
tity and the variance of this measure. Thus we have 

("Las.(0o)) - (^(flo)) ^ 2 - (^sample)) 

<« 2 (0o)) </c 2 (0„)> ' 

A« cas .(0o)) / rj- rr 

-l— 2 , a \ \ = V C22 \ K Sample);; 

( S f cas -) - S 3 _ 
S3 

_o£2l\ (« 2 (Sample)) 6 2 

' say* (k 2 (0 o )) (« 2 (0o)}' 

As moa S . 



(93) 
(94) 

(95) 



Y (if +4C22 - 4f j) (Sample))- (96) 

The numerical calculations give the values of the required 
coefficients c pq (Appendix A), from which we have 



RcaJ^M* 2 ^)) 
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The cosmic errors of course decrease with the size of the 
sample. What the previous results show is that the er- 
rors are already quite small when the size of the sample 
is 25 deg 2 and would allow, for a perfectly well known 
source distribution, a determination of the cosmological 
parameters at the 10 to 15% level. 

7. Discussion and conclusions 

Using Perturbation theory techniques, we have calculated 
the variance and the skewness of the local filtered con- 
vergence k in an inhomogencous cosmology in the FRW 
framework. The quantity k can be calculated from the 
measurement of the magnification // and/or the distortion 
S. 

The exploration of the dependence of the variance of 
k on the parameters of interest have shown that it is, as 
expected, strongly dependent on O , but also on A. It is of 



course proportional to the magnitude of the power spec- 
trum of the matter fluctuations. These dependences, al- 
though degenerate with each other, are of immediate cos- 
mological interest. Note that these dependences are given 
for a given power spectrum normalization at redshift zero. 
As pointed out in other studies (Kaiser 1996, Jain & Sel- 
jak in preparation) if the power spectrum is normalized 
to the present number density of clusters the degeneracy 
with £1 almost vanishes. More worrying is the fact that 
the variance is also strongly dependent on the redshift of 
the sources z s , making difficult to analyze quantitative re- 
sults in a cosmological context if the sources are not well 
known. 

It is possible to disentangle the dependence of the vari- 
ance on the magnitude of the power spectrum and the 
cosmological parameters by the use of two different popu- 
lations of sources. Thus, the ratio of the variances for these 
two populations could provide constraints on the cosmo- 
logical parameters independently of the power spectrum 
(the dependence on its shape is only weak). Of course the 
problem mentioned previously for the knowledge of the 
population of sources is even more critical since both pop- 
ulations have to be perfectly known. 

More interestingly, the skewness S3 of the local con- 
vergence, expressed as the ratio of the third moment and 
the square of the second, can also provide a way to sepa- 
rate the dependence on the cosmological parameters. We 
indeed found that this quantity does not depend on the 
amplitude of the power spectrum, is weakly dependent on 
its shape and varies almost as the inverse of the density 
parameter, with a slight degeneracy with the cosmological 
constant. Moreover the skewness is a probe of the initial 
Gaussian nature of the density field. If physical mecha- 
nisms such as cosmic strings, topological defects, symme- 
try breaking, etc... generate non Gaussian features in the 
initial density field, they will modify the behavior of S3 
(see Gaztanaga & Mahonen 1996), inducing a strong an- 
gular variation of S3 . The dependence of S3 on the redshift 
of the sources is surprisingly large, stressing the absolute 
necessity of knowing the population of sources used to do 
such measurements. It is however interesting to note that 
the product of the variance and the skewness is almost 
independent of the redshift of sources, so in any case we 
have a robust estimator of the product P(k) fl 5 . It might 
reveal to be an interesting quantity to consider. 

An important concern was the contribution of mul- 
tiple lenses configurations in reducing the significance of 
these statistical estimators. We have shown that the two 
main quadratic couplings between two lenses arc inef- 
fective, with contributions of a few percent, and do not 
strongly depend on the cosmology. Thus, whatever the 
cosmological model, the cosmological signal contained in 
the variance and the skewness could be slightly weakened, 
but is not washed out by these additional couplings. 

We have not considered higher order moments, basi- 
cally because we think that they are not going to improve 
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the situation, i.e. to raise the degeneracy between the cos- 
mological parameters we have with the second and third 
moments. A rough calculation (taking into account only 
the dominant contribution) shows interestingly that the S4 
to s| ratio is expected to be close to 2 for any cosmological 
model (l^Oj A and the shape of the power spectrum). In 
principle it could be used as a test for the gravitational 
instability scenario with Gaussian initial conditions. We 
are however far from a reliable measurement of all these 
quantities. 

A concern that has not been addressed in this paper 
is the validity domain of these results with respect to the 
nonlinear corrections due to the dynamical evolution of 
the cosmic fields. These corrections are expected to in- 
tervene at small angular scale where the physical scale 
of the lenses is expected to enter the nonlinear regime. 
The question is naturally, at which angular scale? The 
answer is not straightforward, and actually the situation 
is slightly different for the second moment and for the 
skewness. Indeed the second moment is always a measure 
of the projected power spectrum whether or not it is in 
the linear regime. This point was clearly made in a recent 
work (Jain and Seljak, in preparation) where the rms po- 
larization is calculated including the nonlinear evolution of 
the power spectrum using prescriptions as the ones devel- 
oped by Hamilton et al. (1991), Jain, Mo & White (1995) 
and Peacock & Dodds (1996). For s 3 there are no such 
transformations that would take into account the nonlin- 
ear effects. But one should have in mind that this parame- 
ter has been observed to be very robust against nonlinear 
corrections (this is not necessarily the case for the prod- 
uct of the variance and S3). It has been amply checked 
in numerical simulations for the 3D filtered density field 
(see the remarkable results of Baugh, Gaztafiaga & Ef- 
stathiou, 1995). The situation for the projected density 
field has not been investigated in as many details, and 
extended numerical work would be necessary for that re- 
spect. In particular it would be interesting to extend the 
work of Colombi et al. (1996), who propose a prescription 
to correct the high order moments of the local 3D filtered 
density PDF, to the projected density. Another arduous 
line of investigation would be to take into account in the 
analytic calculations the so called "loop corrections" as 
it has been pioneered by Scoccimarro & Frieman (1996) 
for the un-smoothed density field. Such calculations could 
quantify the errors introduced by the nonlinear effects at 
small angular scale. 

In the last section we have briefly explored the obser- 
vational requirements for completing a large observational 
program. We have in particular tried to estimate the cos- 
mic noise that would affect a catalogue covering an area of 
25 deg 2 . This calculation is at most indicative and should 
not be considered too realistic. For instance, in this anal- 
ysis we have made a fairly crude hypothesis for the noise. 
In particular we have assumed that the noise due to the 
intrinsic ellipticities of the galaxies is local. If this is true 



for the local distortion measurement, this is not true for 
the convergence if it has been obtained from the recon- 
struction scheme proposed by Kaiser since it involves the 
gradient of n. In this analysis we have also not taken into 
account the fluctuations in the number density of sources, 
that could be significant at the angular scale we are inter- 
ested in. 

On the other hand we have assumed that the moments 
were calculated directly, i.e. with the average of the proper 
power of the different measured convergences. This may 
not be the most robust way to do it. Indeed, taking ad- 
vantage of the fact that we expect the Probability Dis- 
tribution Function of the local convergence to be close 
to a Gaussian distribution it is possible to estimate the 
variance and the skewness simply from the shape of the 
PDF around its maximum. More precisely one can use 
the Edgeworth expansion (see Juszkiewicz et al. 1995 and 
Bcrnardeau & Kofman 1995) to look for the best fitting 
values of the low order moments without having to actu- 
ally compute them. With such a method the results should 
be less sensitive to the cosmic noise, because less sensitive 
to the rare events in the tails. 

We leave for a coming paper the difficult task of defin- 
ing what we think could be an optimal strategy for observ- 
ing the shear induced by the large-scale structures and ob- 
taining constraints of cosmological interest out of it. But 
we stress in any case that such a project demands a lens- 
ing survey at 25 square degrees scale with a very high 
image quality. Such a lensing survey will be used for the 
statistical work described in this paper, but also to probe 
the angular power spectrum at scales from the Mpc to few 
hundred Mpcs. This program is accessible with the UH8K 
camera or with the future 16Kxl6K (MEG AC AM) cam- 
era at the Canada-France-Hawaii Telescope, provided op- 
timal algorithms for measuring very weak shear are used. 
The pixel to pixel autocorrelation function (ACF) pro- 
posed by Van Waerbeke et al. (1996) which consists in 
analyzing the cllipticity of the ACF induced by the gravi- 
tational shear seems a promising approach because it does 
not require any more to define object centroids and to 
compute shape parameters for each individual faint galaxy 
from noisy images. The shape of the sample also impacts 
on the signal to noise: reconstruction method a la Kaiser 
is expected to introduce finite size errors along the edges 
of the catalogue, so it would be interesting to have a cata- 
logue as compact as possible. On the other hand an elon- 
gated catalogue has less cosmic variance. 

The drawback of this project is that it requires the 
knowledge of the redshift distribution of the sources, es- 
pecially for faint objects (J5 > 25). A spectroscopic sur- 
vey to get the redshift of these objects is unrealistic with 
present-day ground based telescopes; fortunately many 
other methods are possible, such as the photometric red- 
shifts, the lensing inversion technique, and the depletion 
method (see Mellier 1996 for a review). None of them have 
proved their ability to get secure redshifts, but they are 
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potentially very attractive because they are not time con- 
suming. 
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Appendix A: Expressions of joint cumulants 

General Expressions 

The statistical quantities described in the main text are 
expected to be affected by a cosmic variance due to 
the finite size of the sample. The estimation of such 
an effect is determined by joint moments between local 
smoothed convergences taken at two different locations. 
Indeed these moments provide information on how ob- 
served local smoothed convergences are expected to be 
correlated. The exact derivation of the joint moments of 
interest would be a difficult task, that requires an explicit 
calculation of the higher order expression (up to the sixth 
order!) of the local convergence. As we do not require a 
very accurate calculation but rather a realistic estimation 
of these quantities we will make some simplifications in 
this derivation. First of all we will assume that the high 
order terms of the local convergence are dominated by 
the projection along the line-of-sight of the corresponding 
high order terms of the local density in the dominant lens 
plane. This is what has been found in Section 5 for the 
skewness, and we expect it to be also true for the cumu- 
lants we are interested in. Moreover to do the calculations 
we will assume that the size of the sample in which the 
local convergence is available is significantly larger that 
the smoothing angular scale. This will make the deriva- 
tion of the joint cumulants more simple and accessible 
to Perturbation Theory. Similar results were obtained in 
a slight different context (Bernardeau 1996a) for which 
the derivation of joint cumulants is given for 3D top-hat 
smoothed density contrasts. It is shown that when the 
distance between the smoothing cells is large compared to 
the smoothing scale, the joint cumulants can be written, 

(SfSl)^ C p , q (6!) p - 1 (S^y- 1 (6,6,), (101) 

at the leading order in Perturbation Theory, and where 
the C pq coefficients are finite numbers (they are constant 
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for a power law spectrum). A crucial property for these 
coefficients is that we have, 



where 



C M — Cp. \ Cq.l 



(102) 



so that the derivation of these coefficients reduces to the 
calculations of a single index series of numbers. They have 
been given in the paper mentioned above for the 3D top- 
hat smoothing. 

These results cannot be used directly to estimate the 
joint moments we are interested in, but form however the 
basis of the calculations. Indeed if we assume, as men- 
tioned in the introduction of this appendix, that the high 
order expressions of the local convergence are dominated 
by the projection of the local density of the lenses at the 
same order, then the expressions of the joint moments are 
given by the integration over the line-of-sight of the re- 
lation (101) where the coefficients C Pi q correspond to the 
top-hat smoothed joint cumulants for the 2D dynamics. 
This is a situation which is similar to the situation en- 
countered in the derivation of the expression of the high 
order moments in angular survey. 

Using the calculation developed in Bernardeau (1996a) 
and applying it to the 2D dynamics we get, for a power 
law spectrum of index n, 



2D _ 24 (n + 2) 
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Then, as a consequence, the local joint cumulants fol- 
low a similar hierarchy 

<««i) p «« 2 )*) c = (105) 

^,<«(€i) 2 ) p " 1 («(€ 2 ) 2 ) a " 1 <«(€i)««2)). 

where the coefficients c Pj(? are given by ratios of integrals 
along the line-of-sight. We give here the expressions of 
these integrals for an Einstein-de Sitter Universe, and we 
will limit the explicit numerical calculations of these ex- 
pressions for this case, 
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Jo is the Bessel function of the first kind and 812 is the 
angle between the two directions of the smoothing cells 1 
and 2. 

The final step is to take the geometrical averages of 
these joint moments in the sample. It is expected to be 
dominated by cells being at the distance of the order of 
the size of the sample, # S ampic- Then the averaged moments 
follow the same hierarchy as in (105) but for which the 
function is replaced by 



^2 D {V) 



-F 

Jo 



k ± dk ± P(k±) 



Wi D [k±V9 lt 2] 



(110) 



where W2D is introduced in (45). This result is valid if the 
sample has a spherical shape, although it is by no means 
a crucial hypothesis. 

In the latter case the integral Q 2 is actually the ex- 
pression of the variance of the convergence in the whole 
sample. All the considered joint cumulants are also found 
to be proportional to this variance, and thus the cosmic 
uncertainties on the results will be all the more impor- 
tant that this integral is large. The results will therefore 
depend very strongly on the adopted shape of the power 
spectrum. In the following subsection we give the numer- 
ical results for the APM measured power spectrum (5). 
Compared to a CDM spectrum it contains more power at 
large scale and is thus less favorable for the estimation of 
the cosmic errors, but it is probably more realistic. 

Numerical Results 

In the following we give the results of the previous inte- 
grals for the power spectrum given in (5). The index at 
the 30' smoothing scale is found to be 

nw-1.45, (111) 

and then the coefficients of interest are found to be 



C2,i = 36.7 
c 2 .2 = 1400 
c 3 ,i = 3200 
c 3 , 2 = 125,300 
c 3 , 3 = 11,400,000 



(112) 
(113) 
(114) 
(115) 
(116) 



for a sample size of 25 deg 2 , corresponding to a sample 
size of Sample w 3 deg. These results are calculated with 
the distribution of lenses (2) in which z s w 1. 



(108) Appendix B: Evaluation of the Cosmic Errors 

In this appendix we present the derivation of the errors, 
due to the cosmic variance, i.e. the fact that measurements 
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made in a finite sample are affected by systematic and 
random errors. We will naturally focus our presentation 
on the second and third moments of the local smoothed 
convergence. 

To present the calculations it is necessary to introduce 
two different averages. One is the geometrical average, the 
mean value of a given observed quantity at different loca- 
tions, and the second is the ensemble average, which corre- 
sponds to the expectation value of quantities that depend 
on cosmic variables (local densities..) or random quantities 
associated with intrinsic properties of the galaxies such as 
their intrinsic cllipticities. In principle the two averages of 
a same quantity coincide in a perfect case, that is if the 
available data set was infinite, but in practice this is ob- 
viously not the case! Actually, the geometrical averages, 
considered as estimations of the corresponding ensemble 
averages, are also random variables that are expected to 
vary from one sample to another. In the following we es- 
timate the properties of these geometrical averages, that 
is their cosmic expectation values and their variances, as- 
suming that for a sufficiently large sample they obey a 
Gaussian statistics. We follow here ideas that were devel- 
oped by Colombi, Bouchct & Schacffer (1995) and more 
precisely by Szapudi & Colombi (1996). 

Notations 

In the following we denote X° the (connected) geometrical 
average of a quantity X (and as usual, (X) its ensem- 
ble average) . Thus the observed moments of the measured 
convergence are 



1 N " 
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(117) 
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where the summations are made over the iV c different lo- 
cations where the smoothed convergences are supposed to 
have been determined. We actually assume that the local 
convergences have been measured in a compact sample, 
of a size significantly larger than the smoothing length. 
In the geometrical averages it is also not excluded that 
the smoothing areas overlap, so that N c can actually be 
arbitrarily large. Of course the resulting averages cannot 
be arbitrarily accurate because the measured Ki are cor- 
related variables, and will be all the more correlated that 
they are measured in close directions. 

Modelisation 

In the following we not only estimate the errors due to the 
cosmic variance but also the effects of the intrinsic noise 
in the measurements due to the use of a limited number of 
tracers for the shear and of their intrinsic imperfections. 



So assume that the measured convergence in the direction 
£i is given by, 



«meas.(£i) = + e (£i)> 



(119) 



where is the true cosmological value of the local 

convergence in the direction £ 4 and e(£J is the error made 
in this same direction due to the intrinsic cllipticities of 
the galaxies. Assuming the ellipticities of the galaxies are 
independent from one another, and using results obtained 
for their intrinsic ellipticities (see for instance Miralda- 
Escude 1991b) we get, 
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where n g is the number of galaxies in a given smoothing 
angular cell. Moreover two variables e(£j) and e(£j) are 
independent if the corresponding cells do not overlap, and 
they are all assumed to be independent of the true values 
of the local convergences. This is probably a rather simple 
and naive modelisation but we leave more accurate numer- 
ical studies of a more realistic modelisation for a further 
paper. 

As the variance of e(£j) is independent of i we will 
denote, e their common RMS fluctuation, 



m 2 ) = e 2 



(121) 



and from the expected density of galaxies in a (30') 2 field 
(about 10 4 for the limit magnitudes B = 27 or I = 26 
currently used in ultra-deep imaging; Tyson 1988, Smail 
et al. 1995) we can estimate that 
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We then define the random variable e p that describes 
the departure between an observed geometrical average 
and the corresponding ensemble average, 



«(0o) = ei 
and 



«p(0 o ) = (k p (0 )) c (1 + e p ) 



(123) 



(124) 



for p > 1 . In the following we will assume that the cosmic 
random variables e p have a small variance compared to the 
variance at the smoothing scale. Therefore the variables 
e p can be considered to follow a Gaussian statistics and 
we will do the calculation at their dominant order. The 
knowledge of the cosmic errors in any of the considered 
measurements can then be derived from the values of the 
variances and cross-correlations between these variables. 
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Statistical properties of the variables e p 
To obtain the statistical properties of e p one can sim 
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ply consider the ensemble average of the k p (9 ) , of their 
squares and products. From the assumed properties of the 
variables e(£J and k(^) we have, 



From these results we have entirely define the statistical 
properties of the variables ei, ei and e 3 . Simple identifi- 
cations lead to 
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where we have identified the geometrical average of 
( K (£i) K (£j)) w ith the variance of the local convergence at 
the scale of the sample, Sample- At this scale we assume 
that the ensemble average of e(^) is negligible compared 
to the one of k: it decreases like Poisson noise, that is more 
rapidly than the variance of the convergence. 

The ensemble averages for the second moment read, 
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in terms of the random variables e p and expanded up to 
where the calculations have been limited to the quadratic the quadratic order. This is a trivial calculation for the 
terms in e and (K 2 (6> samp i )) 1/2 . The calculations make variance, and for the measured s 3 we have 
also intervene the geometrical averages of the joint mo- 
ments (K 2 (£i) K 2 (£j)) which are evaluated in the previous 
appendix. 

For the third moment we have, 



It is interesting to note that e enters only in the expecta- 
tion value of e 2 . 

As a result we can calculate both the expectation value 
of a measured quantity and the variance of this measure. 
To do so the quantities we arc interested in arc expressed 
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s§(l + 2e 3 -4e2 + e 2 -8e 2 e 3 



10e 2 



.(1)14) 



We have also to consider the cross-correlations between 
those moments which give, 



Of course this expansion is correct only if the sample size 
is larger than the smoothing scale. 

Results for the variance and for s 3 

Then when we apply the previous rules on the variables 
e p and get 



(Kmeas.(^o) ^meas.^o) ) = 
C 2 ,l (k 2 (#o)) (^sample)) , 

(K mcas .(9 ) c «;Las.(^o) c ) = 

C 3 4(k 2 (^o)) 2 (k 2 (Sample)) 



(«Las.(^0)) - (k 2 (6 )) 



^ Av^ (^sample ) } 



(«meas.(#o) K m cas.(#o) ) = 
(« 3 (0 O )) 2 - <K 3 (0 O )> <K 2 (^ample)) + 
+C 3 ,2 (k 2 (# )) 3 (^sample)) - 
-3 C2.1 <K 2 (#o)) 2 (« 2 (Sample)) +£ 2 (k 3 (# )). 



(131) 
(132) 
(133) 



s -)-«3 



(«taeas.(^0)> 
( 



^sample)} i 



„mcas. 
S 3 



2-3^1 

«3 



S3 

(^K (^sample 



)) 



<« 2 (^o)) 
Asf as - 

mcas. 
S 3 



- 2 



Wo))' 



(145) 
(146) 

(147) 
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^(lf +4C22 - 4 |f) (Sample))- (148) 

One can distinguish two different effects of the cosmic 
variance for the uncertainties of the results. There is first 
a systematic shift proportional to the ratio of the vari- 
ances. And secondly a cosmic scatter is expected to ap- 
pear, which is directly proportional to the RMS of the 
fluctuations of the convergence at the sample scale. It is 
thus crucial to have a catalogue as large as possible! 

Using the numerical results of the previous appendix 
we have 



(^.w)-^(^ = _ 013; 



(k 2 (0o)> 



A< cas .(flo) 

(Kmeas.^o)) 



<sr as -) - ^ 



«3 



0.14; (150) 
-0.07; (151) 



A „mcas. 

5 = 0.04. (152) 



The cosmic errors of course decreases with the size of the 
sample. What the previous results show is that the er- 
rors are already quite small when the size of the sample is 
25 deg 2 and would allow, for a perfectly well known source 
distribution, a determination of the cosmological param- 
eters at the 5 to 10% level. Note that S3 appears to be 
less sensitive to the cosmic variance because it is a ratio 
of moments. 



This article was processed by the author using Springer- Verlag 
WTtfX. A&A style file L-AA version 3. 



